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Developing an expression of resonant x-ray scattering (RXS) amplitude which is convenient 
for investigating the contributions from the higher rank tensor on the basis of a localized 
electron picture, we analyze the RXS spectra from the magnetic phases of Ho near the M^^^ 
absorption edges. At the M5 edge in the uniform helical phase, the calculated spectra of the 
absorption coefficient, the RXS intensities at the first and second satellite spots capture the 
properties the experimental data possess, such as the spectral shapes and the peak positions. 
This demonstrates the plausibility of the adoption of the localized picture in this material 
and the effectiveness of the spectral shape analysis. The latter point is markedly valuable 
since the azimuthal angle dependence, which is one of the most useful informations RXS can 
provides, is lacking in the experimental conditions. Then, by focusing on the temperature 
dependence of the spectral shape at the second satellite spot, we expect that the spectrum 
is the contribution of the pure rank two profile in the uniform helical and the conical phases 
while that is dominated by the rank one profile in the intermediate temperature phase, so- 
called spin slip phase. The change of the spectral shape as a function of temperature indicates 
a direct evidence of the change of magnetic structures undergoing. Furthermore, we predict 
that the intensity, which is the same order observed at the second satellite spot, is expected 
at the fourth satellite spot from the conical phase in the electric dipolar transition. 

KEYWORDS: resonant x-ray scattering, Ho metal, amplitude formula, satellite spot signal, 
helical magnetic structure 



1. Introduction 

With the advent of the so-called third generation synchrotron sources, resonant x-ray 
scattering (RXS) has opened a new research field of investigating, for example, the multipole 
ordering phases in the / electron systems due to the electric multipolar transitions. RXS is 
a second order optical process that an incident photon with the energy tuned through an 
absorption edge excites the core electron to an unoccupied level, then the excited electron 
decays into the core level accompanied by emission of a photon. Element, electron shell, and 
chemical valence specific natures of the RXS technique have made it possible to become a 
complement to neutron scattering to elucidate the spatial Fourier transform of the multipole 
order such as quadrupole and octupole as well as dipole moments.^'^"^^'^'^'^^'^^'^^ 
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Historically, however, analysis of the magnetic properties on the magnetic materials in 
terms of the x-ray resonant process, so-called resonant magnetic x-ray scattering (RMXS), 
developed in its own right since Blume suggested a possibility of such a usage. ^ The first 
observation of RMXS was made by Namikawa et al. on ferromagnetic nickel.^ Since then, 
further studies followed, investigating, in particular, the / electron systems such as holmium, 
uranium arsenide, and so on.^°'-^^'^^ 

Metallic holmium is one of the most studied materials by using RMXS. The crystal struc- 
ture of holmium is hexagonal close packed structure with lattice constants a = 3.58 A and 
c = 5.62 A at room temperature. On cooling down the temperature, magnetic moments in 
Ho order at Tn = 133 K with a helical magnetic structure, which consists of the ferromagnet- 
ically aligned moments within each basal plane and the orientation of the moments rotates 
in successive basal planes. As the temperature decreases, the wavevector characterizing the 
spatial modulation of the helical pattern reduces and locks in below Tc = 20 K where the 
magnetic structure becomes a conical configuration.^'^ 

Gibbs et al. detected the RMXS signals near the Ho L3 edge in the magnetic satellite 
spots. Polarization analysis distinguished two kinds of spectra with peak positions separated 
about 6 eV. They are attributed to the consequences of the electric dipole (i^l) and quadrupole 
{E2) transitions. The former had intensities up to the second harmonic satellite spot while 
the latter was up to the fourth harmonic satellite spot. Hannon et al. gave the theoretical 
explanation of the experimental results. They derived useful formulae describing the RMXS 
amplitude and, by using them, concluded that relatively large scattering amplitude would 
be expected in the vicinity of the Ho L2.3 and M^^ edges, the former of which explained 
the experiment. By extending the formula derived by Hannon et al., Hill and McMorrow 
investigated some qualitative properties of the RXS signals at the Ho L3 edge in the uniform 
helimagnetic and the conical phases. These works have stimulated the immense activities in 
the field of RMXS both experimentally and theoretically.^'^ 

As for Ho, several experiments followed focusing, for instance, on the temperature de- 
pendence of the exponents, the existence of two-length scales, and so on.^^'^^'-*^^ Theoretical 
analyses based on the mean field model revealed the importance of the interaction having the 
trigonal symmetry. Despite the numerous achievements provided by the RMXS data at 
the Ho L3 edge, it is difficult to extract the quantitative interpretation since the El process 
at the L3 edge is the consequence of the transition between the 2p and Be? states so that the 
relevancy of the 4/ electrons is indirect. The band nature of the hd electrons also increases the 
complication of the analysis. Although the E2 signal, the transition between the 2p and 4/ 
states, presents a direct information of the 4/ electrons, the interference with the El signal 
complicates the quantitative interpretation. 

A possibility has long been suggested that the RMXS measurement at the Ho M4^5 edges 
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by means of soft x-ray beam may supply the useful insights to elucidate the natures of the 
magnetic structures. ""^^ The El processes at the M^^^ edges are brought about by the transition 
between the 3d and 4/ states, which provides a direct information from the 4/ electrons. 
Furthermore, long wavelength of the soft x-ray is suitable for investigating the long period 
magnetic structure of Ho. Recently, Spencer et al. and Ott et al. have succeeded in detecting 
the RMXS spectra of the first and second harmonics satellite spots as well as the absorption 
coefficient spectrum at the Ho M5 edge from the helical magnetic phase. ^^'^^ The energy 
profiles of the spectra of the first and second satellite spots differ significantly from each 
other. They were attributed to the circular and linear dichroic contributions, respectively. In 
this case, it turned out that the absence of the E2 contribution made the analysis rather 
simple. These facts make a sharp contrast observed at the L3 edge where the difference of the 
spectral shapes at the different satellite spots were not so prominent. 

Usually, the azimuthal angle (tp) dependence of the peak intensity of the spectrum offers 
one of the most highlighted outcomes given by the R(M)XS experiment in many cases. In 
holmium, however, the intensity shows no ip dependence in the uniform helical and conical 
phases with the choice of the scattering vector adopted in the experiment. Therefore, we should 
extract implications as many as possible through the spectral analysis. Then, a theoretical 
study may lend support to the analysis of the R(M)XS spectra at the Af4.5 edges. 

In this paper, therefore, we perform a theoretical spectral analysis of RMXS at the Ho M4^5 
edge. Based on a localized electron picture, we calculate the spectra exploiting a useful expres- 
sion of the scattering amplitude whose utility is examined in several / electron systems. 
We investigate three spectra; the RMXS spectra at the first and second satellite spots and 
the absorption spectrum. Their shapes at the M5 edge are distinguishable one another and 
capture the main features presented by the experiment, spectral shape and the shift of the 
peak position at the different rank satellite spots. The ratio of the peak intensity extracted 
from the calculated spectra also is the same order as the experimental one. 

Another result we present is a prediction that the spectral shape varies as a function of 
temperature. For instance, the spectrum at the second satellite spot shows the same energy 
profiles both in the helical and conical phases. On the other hand, the profile becomes com- 
pletely different one in the temperature range between two phases, known so-called as the 
spin slip phase. Such an evolution of the spectral shape may provide the direct evidence that 
the change of the magnetic structure is undergoing. 

Further inferred from our investigation is a possibility of the observation of the RXS 
spectra at the higher-order satellite spots. In the conical phase, the presence of fifth- and 
seventh-order magnetic satellites in the neutron scattering were observed.^ This is attributed 
to the fact that the turn angle of the magnetic moment between successive Ho planes is not 
a constant (see Fig. 1 (b)). Similarly, we show that the RXS intensities are expected at the 
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higher-order satelhte spots in this phase. For example, in the El transition, fourth-, fifth-, 
seventh-, and eighth-order satelhtes are possible to be detected. Our estimation reveals that 
the RXS spectrum at the Ho M5 fourth-order satellite spot from the ir — a' channel possesses 
the same order of magnitude and the same energy profile as that observed by Ott et al. at 
the second satellite spot. 

The El transition at M4,5 edges of rare earth metals give rise to an excitation of electron 
from the 3(i levels to the 4/ levels. Since the associated photon energy is about 1 ~ 2 KeV, 
this process is only useful to investigate the long-range orders with rather long periods in 
comparison with the lattice constant. In this paper, taking up the helical and conical magnetic 
orders in Ho metal, we analyze the RXS spectra in detail 

The organization of this manuscript is as follows. In Sect. 2, we explain the various mag- 
netic structures exhibited by Ho. In Sec. 3, we present a brief summary of the theoretical 
framework of R(M)XS. An expression of the R(M)XS amplitude the present authors had de- 
rived^^'^^ is introduced and is compared with the previous similar results. In 
Sec. 4, the initial and the intermediate states needed to calculate the RXS spectra are defined 
in order to apply the theory to the investigation of the RMXS spectra of the satellite spots in 
the vicinity of the Ho M5 edge. Qualitative features expected from each phases are analyzed 
in Sec. 5. Numerically calculated spectra of the absorption coefficient and RMXS are shown 
and compared with the experiment in Sec. 6. Main conclusions are summarized in Sect. 7. 

2. Magnetic ordering phases of Ho metal 

Before calculating the RXS spectra, we summarize the magnetic structures experimentally 
identified in the Ho metal. The crystal takes an hep structure. Below Tn = 133 K, the magnetic 
moments of the system order in a helical structure. The moments on each Ho layer are confined 
to the basal plane and are coupled ferromagnetically within each layer. The interlayer angle 
of the moment varies a constant angle per layer. We call it as ''^uniform helical phase". A 
schematic explanation of the orientation of the moment on each layer of this phase is given in 
Fig. 1 (a). The number of Ho layers N within one helical pitch is connected to the modulation 
vector q = (0, 0, 2/Ar) in units of 27r/c. The N increases on cooling because the helical structure 
is distorted, or more precisely, the magnitude of q reduces from q = (0, 0, 0.271) corresponding 
to AT ~ 7 just below Tn to q = (0, 0, |) corresponding to AT = 12 at Tq = 20 K. Below TtextrmC, 
q is locked in to (0,0, g) corresponding to A" = 12 with the emergence of a ferromagnetic 
component along the c-axis, called as " conical phase" . A complication is that the orientation 
of the moment within a basal plane does not alters constantly.^' ^ 

The magnetic easy axes seem to direct to the centers of six triangles formed by neighboring 
atoms. The actual structure in the conical phase is that the twelve Ho layers form six pairs of 
two layers with the moments directing nearly to one of the centers of six triangles formed by 
nearest neighbor atoms as shown in Fig. 1 (b). The bunching angle between the pairs across 
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the same easy axis is defined as 7 and evaluated about 5.8° in the T ^ hmit. 

The transition between the uniform helical and the conical phases is not a straightforward 
depending on the conditions such as the strain present in the sample, the presence of the 
chemical impurities, and so on.-*^^ When the temperature is above Tq but well below Tn, the 
c component of the moment disappears but the bunching remains. At the same time, some 
of the pairs lose their partner and the remaining moment directs to one of the centers of 
six triangles formed by nearest neighbor atoms. Gibbs et al. called this phase as the ''spin 
slip phase". For example, when one pair experiences a spin slip, the unit cell of the system 
becomes five bunched pairs plus a single component with = 11 (See Fig. 1 (c)). When 
the temperature raises, the value of 7 increases resulting in the uniform helical phase. In this 
paper, we show that the spectral shapes of the RXS spectra changes drastically when the 
magnetic structure changes from the uniform helical one to the others. 



(a) uniform helical phase (N=8) (b) conical phase 
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(c) spin slip phase (N=11) 
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Fig. 1. Orientation of the magnetic moment projected onto the ah plane: (a) uniform helical phase 
with N = 8, (b) conical phase, and (c) spin slip phase with N = 11. Direction of arrow associated 
with number n = 1, • • • ,N represents the direction of the magnetic moment in the n-th Ho layer. 
The dotted lines denote six equivalent easy axes. 



3. An extension of the RXS formula of Hannon et al. 

A resonant process is described as following: an incident photon with the energy huj, wave 
vector k, and polarization e is scattered into the state with the energy fiio, wave vector k', 
and polarization e' through the process that the core electron is excited into an unoccupied 
level leaving a core hole and then recombines with the core hole by emitting the photon. 
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In a localized electron picture, the scattering amplitude is well approximated by a sum of 
the contributions from each site. Therefore, the RXS amplitude /Ei(e, e', k, k'; a;) in the El 
transition is given by 

/ei = ^5;^e-«-'-^M,(e,e';a;), (3.1) 



with 

a'a A 

within the second order of the photon-electron interaction in the perturbation theory. Here 
the scattering vector is defined as G = k' — k, and the number of sites is represented as Nq. 
The Mj{e, e';uj) is the scattering amplitude of single site j with the position vector rj, where 
the initial state is represented as |0j) with the energy Eq and the intermediate state is given by 
I A) with the energy E\. The lifetime broadening width of the core hole is denoted as F. The 
dipole operators XaiJYs are defined as xi{j) = x, X2{j) = y, and x-i(j) = z in the coordinate 
system fixed to the crystal axes with the origin located at the site j. The RXS intensity is 
proportional to the square of the absolute value of the scattering amplitude. 

The evaluation of eq. (3.2) is not an easy task, because the intermediate states are to be 
summed up. We consider the situation that the initial state is described within a multiplet 
of certain total angular momentum J, and the degeneracy is lifted by the crystal electric 
field (CEF) and the inter-atomic interaction. This situation is familiar in many f electron 
systems. On the other hand, we neglect the CEF and the inter-atomic interaction on the 
intermediate states. This may be justified because the energy dependence of the spectra is 
mainly determined by the multiplets which energy is mulch larger than the CEF and the inter- 
atomic interaction. With assumption, the intermediate states possess the rotational symmetry, 
and thereby the summation over these states are rather simply carried out. We could express 
the scattering amplitude at site j in a simple form as 

Mj = ag(a;)(e'-e)-ia£(a;)(e'xe)-(0,|J|0,) 

+ ag(a;)^PP(6,6')(0,|i.|0,). (3.3) 
i/=i 

Here (a;) represents the energy profile of rank n operator, which depends only on the initial 
state through the magnitude of the angular momentum J. The geometrical factors associated 
with the rank two quantities are represented as PiP{e,e'). The components of the (rank two) 
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quadrupole operator are defined as 

d3,2_,2 = i[3j2- J(J + 1)], V = 2 

dyz = ^iJyl + jJy), 1^ = 3 . (3-4) 

Ozx — ^^{Jz'Jx ~\~ •Jx'Jz)^ V = A 

Oxy — ^^(^*^x'^y ~t~ *^y^x^^ = 5 

The detailed derivation of the above expression including the definitions of quantities (a;) 
andPi^^(€,e') are found in our previous papers,^^ in which applications to several / electron 
systems and an extension to the E2 process are explored too. 

The formulae similar to eq. (3.3) are found in the previous literatures.^^' ^^'^^'^^ One of 
the most notable formulae is the one derived by Hannon et al.,^^ which is expressed for the 
magnetic ordered system as 

M, = (e' • e)/^(a;) - i{e' x e) • ni,f[{u;) + (e' • m,)(e • m,)/^(o;), (3.5) 

where nij is the unit vector in the direction of the magnetic moment at site j. The energy 
profiles are represented by /'j(cl')'s describing the contributions from the order n of the mag- 
netic moment. Although both arc similar in a symmetrical point of view, our expression is 
more accurate and superior than that of Hannon et al. in the following reason. 

First, higher rank multipole is naturally treated by our expression, while that of Hannon 
et al. was expressed in terms of the expectation values of the dipole moment operator. To 
make the comparison easy, we rewrite the last term of eq. (3.3) as 



(2)/ \^ 
"El 1^)2 



i((6 • J)(6' • J) + (e' • J)(6 . J)) - i^:l±}l(^e' ■ e] 



(3.6) 



In comparison with Hannon et aVs, the expectation value of the product of the operators 
is generally different from the product of the expectation values of the operators. Our ex- 
pression is applicable to other multipole ordered phases without the local magnetic moment. 
The numerical amount of the difference between eqs. (3.3) and (3.5) depends on the system 
considered in the magnetic ordered phase. In the present case of holmium, the difference is rel- 
evant to the analysis of the second harmonics satellite spot, which brings about a quantitative 
outcome when we use eqs. (3.3) and (3.5) as a fitting function. 

Second, our expression gives the spectral shape of the RXS intensity as a function of the 
incident photon energy, without relying on the so-called fast collision (FC) approximation. 
Recent development of the high resolution in energy enables the detection of the reliable 
spectral profiles as a function of the incident photon energy, which increases the value of our 
treatment including the energy dependence. 
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Before going to the detailed analysis of RXS spectra in Ho, we briefly explain how to 
evaluate eq. (3.3) based on a localized electron picture. 
4.1 Evaluation o/(0j|J|0j) and {Oj\zi^\Oj) 

The holmium in the solid behaves, to a good approximation, as trivalent ion with the 
(4/)^^^ configuration equivalent to in the hole picture. We expect from Hund's rule the 
states of ^Ig as the ground multiplet, which is denoted as | J = 8, J^). The degeneracy of this 
multiplet would be lifted by the intersite interaction, resulting in the magnetically ordered 
states. Therefore, the magnetic wavcfunctions at each site are expanded in terms of |J = 
8, J^)'s. It is known that the local magnetic moment is as large as about 10;Ub in the conical 
phase. This value is equivalent to the maximum one expected from the saturation value of 
^Ig. It is also known that the local magnetic moment is about 70 % of the saturated value 
around T/T^ ~ 0.6 in the uniform helical phase. Accordingly, we assume that the magnetic 
moment at each Ho site is saturated, since the spectral shape of the energy profile a^^{uj) 
is not sensitive to the value. Therefore, the wavefunction is given by \ J = 8,Jz = 8) at each 
site in the local coordinate frame where the z-axis is pointing to the direction of the local 
magnetic moment. 

Let (^?„,/3, 0) be the Euler angle connecting the local coordinate frame to the crystal 
fixed frame for the n-th Ho layer, and be the corresponding magnetic state. Then, the 
expectation values of the dipole and the quadrupole operators defined by the crystal fixed 
frame are evaluated as 
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(4.2) 



These expressions are used for (Oj|J|Oj) and (OjIzi/lOj) in eq. (3.3). 
4-2 Evaluation of a^^i{uj) 

Energy profiles a^^-^{uj) {i = 0,1,2) depend on the initial state only through J of the 
ground multiplet but little on the details of the magnetic orders. They also depend on the 
intermediate states. For the ground multiplet, we consider the intra-atomic Coulomb and the 
spin orbit (SO) interactions, where the parameters needed to evaluate the above interactions 
such as the Slater integrals for the Coulomb interaction and the SO coupling constant are 
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calculated within the Hartree-Fock (HF) approximation J' '^'^ Representing the Hamiltonian 
with the bases spanned by the configuration, wc diagonalize the Hamiltonian matrix. We 
obtain the degenerated lowest seventeen states corresponding to J = 8. In the intermediate 
states, the electronic configuration becomes {3d)^ (Af)^^ equivalent to d^f^. We repeat the 
similar procedure as have done to prepare the ground multiplet. The Hamiltonian describing 
the intermediate states takes full account of the intra-atomic Coulomb and the SO interactions 
among the (3d)^(4/)^^ configuration. By diagonalizing numerically the Hamiltonian matrix 
represented with the bases spanned by the d^f^ configuration. We obtain the set of the 
intermediate states. These states are sufficient to evaluate the energy profiles a^^^{u}). 

5. Analysis of RXS for Ho 

Since all the Ho atoms within the same layer give a same amount of contribution, the 
summation over site j is replaced by that over n. The amplitude for the m-th satellite spot is 
given by 

TV 



/ei oc (e • e)aEi (w) ^ e 



n=l 

TV 

..(1) 



.n=l 

N 



,n=l 



(5.1) 



where the wavevector of the m-th satellite spot is defined by 



2 

Gm = (0, 0, Tm) , Tm = m— (5.2) 

which is measured in units of The position vector of the n-th Ho layer is represented as 
r„ = (0,0,§n). 

Now we describe 6n in various magnetic ordered phase. In the conical phase, it is given by 

J 27r;^ + ^0 — 7 for n = odd 



27r^^yi + ^0 + 7 for n = even 



(5.3) 

N 

where is a constant and taken to be zero. The number of the holmium layers within a 
helical pitch is A?^ = 12. In the uniform helical phase, we could put N as an arbitrary integer, 
00 as an arbitrary angle, and 7 = in eq. (5.3). Thereby we have 6n for an arbitrary integer 
n, 

277,-1 

On = + <^o- (5.4) 

In the spin slip phase, we take a case of = 11 system as an example, since the distribution 
of On's show a wide diversity. The unit cell consists of the five pairs of bunched doublet 
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and a single layer as schematically shown in Fig. 1 (c). Within the one helical pitch (for 
n = 1, 2, • • • , A''), 6nS are expressed as 



2vr^ + (/)o-7 forn = 1,3,--- ,Af-2 
2vr^ + (po + j forn = 2, 4, • • • , - 1 



(5.5) 



27r 



for n = iV 



N+l 

with (^0 = f • Note that eq. (5.5) is valid only forn = 1, • • • ,N, which is enough for the present 
analysis. 

With the use of eqs. (4.1) and (4.2) together with eqs. (5.3), (5.4), and (5.5), the calculation 
of eq. (5.1) results in evaluating the summations of the types 
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where £ is an arbitrary integer. Note that eq. (5.6) has no contribution when ^ is a odd integer 
in the uniform helical phase, since the cosine factor becomes cos (^) with 7 = In the case 
of the conical phase, eq. (5.6) is applicable for m' > 1 and when m' = 0, the summation is 
trivial as E^=i e-^^-""" = NSmfi- 
5.1 Uniform helical phase 

In the uniform helical phase, the observed value of N is larger than about seven. Then, 
only the term proportional to S^_^, m with £ = in eq. (5.6) is relevant for the description 
of the El process in this phase. Practically, eq. (5.6) is rewritten as 
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By substituting eq. (5.7) into eq. (5.1), and fixing /3 = 7 = |, we obtain the RXS amplitude 
at the first satellite spot as 
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(5.8) 



Note that the contributions from the quadrupole operator disappear in the uniform helical 
phase at the first satellite spot. Similarly, the amplitude at the second satellite spot is given 
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by, 

« f '^-^<y-" eM^-w)[pf) - iPP]4l(.), (5.9) 

where arguments are omitted in the geometrical factor ■ As aheady known/°'^^'^^ the 
spectrum at the m-th sateUite spot consists of pure rank m profile. The final forms of the 
scattering intensities are derived by substituting the geometrical factors into eqs. (5.8) and 
(5.9). We relegate both the detail of the derivation and the results to Appendix . Here, we 
mention an important feature the results possess that the intensities exhibit no ip dependence 
at both satellite spots with the choice of the scattering vector as the form in eq. (5.2). 
5.2 Conical phase 

In the conical phase, N and 00 are fixed to be twelve and zero, respectively. A notable 
difference compared with the uniform helical phase is that the terms with £ = ±1 becomes 
to be relevant in eq. (5.6) in the conical phase. Let us state the result with £ = and ±1, in 
turn. 

First, the relevant terms in eq. (5.6) with £ = are summarized as follows. At the first 
satellite spot (m = 1), the amplitude is written as 

/ei oc 6e~*T^ cos ^7 - J sin /3 

-i{{e' X e)^ - i{e' x e)y}a^^l{u) 



(5.10) 



The amplitude is a mixture of a^l (lo) and a^l {lo) , while the one is pure {uj) in the uniform 
helical phase. This indicates that there is a possibility to observe the variation of the spectral 
shape of the first satellite spot when the underlying magnetic structure changes. However, such 
a possibility is practically denied since |ag^^(a;)| is much larger than |ag-^^(a;)| as is confirmed 
numerically in the next section. 

The amplitude at the second satellite spot (m = 2) is given by, 

/ei oc 6e~'t cos ^27 — sin^ (3 

X ^ J(2J - l)[pf ) - iPf ]ag(a;). (5.11) 

Note that both amplitudes do not include a contribution from {9n\Jz\Gn) despite the fact that 
the conical phase is characterized by the non-zero value of it. 

Next, we proceed to summarize the result expected from the terms with ^ = ±1 in eq. 
(5.6). This raises a possibility that signals from the higher-order satellites may be detected. 
The possibility is linked to the distortion of the direction of the magnetic moment within 
the basal plane. Koehler et al. have reported that such higher-order magnetic satellites as 
fifth- and seventh-order spots were present in the neutron scattering.^ In the El process of 
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R(M)XS, we see the possibility that the intensities from the fourth-, eighth-order in addition 
to fifth- and seventh-order satehites are detectable. The possibility of detecting the former 
two satellite enable RXS to be a complement to neutron scattering. 

The fifth and seventh satellite spots have contributions from rank one and two profiles. 
The intensities are proportional to 

,2 



OC cos 



(5.12) 



where m = ±5 and ±7. The upper and lower signs in the equation are taken in accordance 
with that of m. The intensities at the fourth and eighth satellite spots have contribution from 
rank two profile, being proportional to 

OC cos^ (27 + I) |Pf ± ^Pf (5.13) 

where m = ±4 and ±8. The upper and lower signs in the equation are taken in accordance 
with that of m. The detectabilities of these intensities at the higher order satellite spots are de- 
termined by the cosine factors and the geometrical factors appeared in the above expressions. 
If we assume 7 = 5.8°, the former factor gives 0.08, 0.19, 0.48, and 0.64 with m = ±4, ±5, ±7, 
and ±8, respectively. Near Ho M5 edge, the Bragg angle 6'("^) = 33.1°, 43.1°, and 73.0° with 
respect to m = ±4, ±5, and ±7, respectively. With the substitution of these parameter values, 
we can conclude that the intensity from the a^l (w) contribution at the Ho M5 fourth satellite 
spot in the tt — a' (or cr — tt') channel is large enough to be detected. Actually, the intensity 
is expected as nearly the same as that observed at the second satellite spot by Ott et al.^^ 
Finally, we end this subsection with a comment that the intensities present no azimuthal angle 
dependence with the choice of as eq. (5.2) in the conical phase as shown in App. . 
5.3 Spin slip phase 

Since the magnetic structure of this phase is filled with variation, we restrict ourself on 
the case of N = 11 system depicted in Fig. 1 (c). The distribution of 0„'s is irregular at the 
spin-slipped layer in this phase. Then, the summation over n similar to the left hand side 
of eq. (5.6) is not results in the right hand side of it. As a consequence, there remains the 
terms for m ^ m'. Two prominent properties are readily anticipated. First, the intensities at 
the first and second satellite spots exhibit ip dependence. Second, the profiles consist of the 
mixture of a^l (a;) and a^l (w) at both spots, which leads to the situation that both spectra 
are dominated by a^l{uj). The latter point at the first satellite spot is equivalent to the 
situation observed in the conical phase, eq. (5.10). It brings about no practical change of the 
spectral shape since a^l{uj) always dominates the entire spectrum. On the other hand, the 
RXS amplitude at the second satellite spot is completely new one, since those in the conical 
and uniform helical phases are made of pure a^l{oj). 
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5.4 Clockwise or anti- clockwise 

In the helical (or eonical) magnetic structure of holmium sample, usually, there exist 
two domains distinguished by the orientation of the winding of the magnetic moment as 
clockwise and anti-clockwise. The latter is defined such that the angle ^„ of the moment 
increases anti-clockwise when the number of Ho layer n increases as shown in Fig. 1, while 
the former the opposite. Since experimentally observed RXS intensities are expected to be the 
domain averaged quantities, we should mention what is anticipated from the clockwise and 
anti-clockwise domains. 

By checking eqs. (5.8) and (5.9), it turns out that the RXS intensity shows no domain 
dependence in the uniform helical phase. The same is true at the second satellite spot in the 
conical phase. The intensity at the first satellite spot, however, shows domain dependence when 
the scattering vector is not the (0, 0, ri)-type. For example, the intensities with G = (2, l,ri) 
may give domain dependence. Fortunately, such G is impossible to achieve at the Ho M5 
absorption edge. When the experiment is possible for this scattering vector at the L3 edge, 
the domain consideration will be necessary interpreting the result. As a result, it is no need 
to worry about the domain consideration when analyzing the conical and the uniform helical 
phases. As for the spin slip phase, the spectrum shows the domain dependence, at least in 
a qualitative sense, regardless of the scattering vector. So we take domain average when 
numerical evaluation is carried out in the spin slip phase in the next section. 

6. Numerical results 

We are now in a position to calculate the RXS spectra in order to compare with the 
experimental result. Before the comparing process goes forward, we first expect that the 
sample used in the measurement is in the uniform helical phase, not in the conical nor spin 
slip phases, because Ott et al. explicitly did not reported such information nor temperature 
the measurement were performed.^^'^^ The data were obtained from the Ho-metal thin film 
with 16 mono-layers (ML). The critical temperature Tn of the 16 ML sample is about eighty 
percent of the bulk one^^ and the first harmonics satellite peak is observed around ri ~ 0.3 
^-1 38 rpj^-g yg^j^g Qf suggests that the temperature is T = O.BSTn inferred from Fig. 4 of 
ref. 38. In bulk Ho, the ratio T/Tn = 0.55 corresponds to T ~ 73 K and n = 0.252 A"^ with 
N 9. On the other hand, if we apply the value of ri to the bulk leading to T ~ 115 K 
with N ~ 7.5J Since the uniform helical phase is transformed into the spin slip phase below 
about 50 K with ri 2± 0.22A in the bulk, both estimations lead us to a conclusion that the 
experiment were performed in the uniform helical phase of the sample as anticipated. Later, 
we will lend a further support that this assumption is sustained when we calculate the second 
harmonics satellite spectrum. We show if the sample was in the spin slip phase, the spectral 
shape might be completely different one. In the following, we assume N = 8, i.e., four unit 
cells are involved in one helical pitch. We emphasize that the choice of TV does not alter the 
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6.1 Absorption coefficients 

We discuss the absorption coefficient A{uj) in the El transition, which is described as 

i A a ' 

It is obtained from eq. (3.3) by the relation 

A{u) (x-Ima^^l{uj), (6.2) 

where ImX represents the imaginary part of X. Although the absorption spectrum at the Ho 
Ms edge was already investigated by others,^' here, we utilize it in three-fold purposes; the 
determination of the origin of the energy, that of the value of F, and the justification of the 
atomic picture. 

The calculated result is shown in Fig. 2. The origin of the energy is adjusted such that the 
maximum point is located at the experimental one, i.e., at 1353.2 eV and is fixed throughout 
this paper. Since the peak structure of the spectrum shows F dependence, several values of 
F are examined. As shown in Fig. 2, we find F ~ 0.6 eV gives a better choice. We do not 
try, however, to obtain the best fit to the data considering the simple treatment based on the 
atomic picture. 

This value of F is consistent with the previous work for the absorption spectra.^ We confirm 
that the small change of the F value around 0.6 eV has no effect on the peak positions. Our 
result reproduces the characteristic multi-peak structure exhibited by the experimental data.^^ 
The total width of the spectrum is governed by the multiplet splitting of the energy of the 
intermediate states. Aside from the peak position around 1348.5 eV in the experiment, we 
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conclude that the atomic treatment gives a reasonable result. Fortunately, we shall confirm 
in the next subsection that this discrepancy causes no practical effect on the investigation of 
the RXS spectra, since the intensities near 1348.5 eV at the first and second satellite spots 
are weak and merely the tail part of the spectra. 
6.2 First and second harmonics satellite 

The RXS spectral intensities at the first and second harmonics satellite spots are propor- 



tional to the dipole 



(1)|2N 
El I ) 



and the quadrupole Profiles, respectively. The analytical 



forms are derived by substituting P = ^ and 7 = into eqs. (A-4) and (A-5) as shown below. 
The intensity at the first satellite spot is proportional to 
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(6.4) 



As seen from the above expressions, the RXS intensities show no dependence. Note 
that spectral shapes have no polarization dependence. In order to compare our results with 
the experimental ones, we consider the spectra in the tt — a' and tt — tt' channels. The ratio 
jW/jC^) is common in both channels with an aid of the relation sin 0^2) = 2sm9^^\ Thus, 
we can fix ^ = and restrict in the tt — a' channel hereafter. The calculated spectra of I^^^ 
and /(2) aj-e shown in the upper and lower panels, respectively, of Fig. 3 together with the 
experimental data.^^ 

As for the \a^l p from the first satellite spot with Gi = (0, 0, ri), the peak position of the 
spectrum appears at the same position of the absorption coefficient, which is coincide with 
the experiment. The entire spectral shape with T = 0.6 eV agrees remarkably well with the 
experimental one including the kink in the low energy tail and the hump in the high energy 
tail. 

As for the la^^^p from the second satellite spot with G2 = (0,0, r2), the maximum height 
intensity is located at the energy 1.8 eV lower than that of the absorption as the same position 
as the experiment. The entire spectral width and the shape of the tail parts show quite similar 
to the experiment, too. On the other hand, however, a prominent discrepancy can be found 
between the theory and the experiment. The peak around 1353 eV found in the the calculated 
result is absent in the experiment. The reason of this discrepancy is still unclear. Prom the 
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Fig. 3. The RXS spectra near the Ho Afs edge at the first (upper) and second (lower) satellite spots. 
The curves for T — 0.4 eV and 0.6 eV are distinguished by the dashed and solid lines, respectively. 
Open circles represent the experimental result. 



with r = 0.6 eV. This large ratio allows us to expect the dominance of a^l (a;) in the entire 
spectral shape when both profiles are mixed. 
6.3 Mi edge 

We carry out the similar calculations at the Ho M4 edge as those at the M5 edge. The 
absorption coefficient is shown in the upper panel of Fig. 4. It shows a single peak structure 
and the peak height of the M4 spectrum is about 0.175 times the one at the M5 edge. These 
observations agree with the experiment.^ The RXS spectra expected from the first and second 
harmonics satellite spots are shown in the lower panel of Fig. 4. The maximum intensity of 
the second satellite spectrum is found at the energy about 0.6 eV higher than that of the first 
satellite spectrum. This tendency is contrary to the case observed at the M5 edge where the 
former is found at the energy about 1.8 eV lower than the latter. The peak intensities at the M4 
edge are about two times weaker than those at the M5 edge. Since the order of the magnitude 
of the spectrum of the first satellite spot at the M4 edge is roughly the same order as that of 
the second satellite spot at the M5 spot, the former may be detectable experimentally. 



calculated value of the ratio I^^^ we deduce a ratio 




(6.5) 



16/24 



J. Phys. Soc. Jpn. 



Full Paper 




Ho edge 



to 



r=0.4eV 

r=0.6 eV 




CD 




Q. 







TO 



-« 0.4 



c 

13 



G,=(0,0,x,),r=0.4eV 
G,=(0,0,x,), r=0.6eV 
G^=(Ofi,z^), r=0.4 eV 
G^={0,0,z^), r=0.6 eV 





\ 



c 



x100 



^ 
1385 



1389 



1393 



(0[eV] 



Fig. 4. Absorption coefBcient and the RXS spectra at the Ho M4 edge. Upper: the absorption coef- 
ficients with F = 0.6 eV and 0.4 eV are shown by the soHd and dashed Hnes, respectively. Lower: 
the RXS spectra at the first and second harmonics satelHte spots are represented by the thick and 
thin hnes, respectively. The values of F = 0.4 eV and 0.6 eV are distinguished by the dashed and 
solid lines, respectively. 

6.4 spin-slip phase 

We investigate the RXS spectra at the satellite spot in the spin slip phase. The qualitative 
evaluation in Sec. 5.3 indicates that they are the mixture of the profiles a^l{u) and a^l{uj). 
The dominance of the former over the latter leads to the expectation that the spectra at 
the first and second satellite spots in this phase are governed by the former. The situation 
at the second satellite spot is of particular interest since the the spectra consist of the pure 
a^l{io) both in the conical and the uniform helical phases. Quantitatively, we calculate the 
RXS spectrum at the second satellite spot in a spin slip phase with iV = 11. The core hole 
lifetime broadening is set to be 0.6 eV. As shown in Sec. 5.4, the spectra exhibit the domain 
dependence, that is, the clockwise or anti-clockwise alternation of the direction of the magnetic 
moment affect the spectra. Thus we take the domain average of the spectra assuming equal 
probability of the two domain. Moreover, the spectrum varies as a function of i/j in this phase. 
The results are summarize in Fig. 5. Although the spectral shape varies with ijj, the entire 
shape including the highest peak position is dominated by the pure |Q;^\^(a;)p as anticipated 
qualitatively. 

This fact brings us to the speculation of the spectral shape evolution. With decreasing 
temperature from Tn down to Tq, the spectral shape at the second satellite spot evolves from 
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1345 1355 
(o[eV] 

Fig. 5. The RXS spectra at the Ho M5 edge with G2 = (0, 0, T2) in the spin shp phase with N = 11. 
The intensities are averaged between the anti-clockwise and clockwise helical phases. The solid 
and broken lines represent the tt — a' and tt — tt' channels, respectively. The azimuthal angles are 
chosen as (a) tp = 0, (b) tp = j, and (c) f/' = f • 

the pure |ag^^(a;)p, through the |Q;^\^(a;)p dominated one, into the original one accompanying 
the shift of peak position from 1353.2 eV to 1351.4 eV and again to 1353.2 eV, respectively, 
large enough to detect experimentally. An experimental observation of the second satellite 
spot spectra in the different phases is desirable. 

7. Concluding remarks 

Stimulated by the RMXS experiments, "^"^^'^^ wc liave investigated the Ho M5 spectra at the 
first and second satellite spots. The analysis have performed by exploiting the useful expression 
of the scattering amplitude derived on the basis of the localized electron picture.^^'^'' Although 
a large enhancement of the intensity is expected at the Ho M edge in the RXS spectrum, the 
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spectrum shows no azimuthal angle dependence in the uniform helical and conical phases with 
the choice of Gm adopted in the experiment. Thus, we have mainly focused on the spectral 
shape analysis. 

A qualitative analysis have revealed that the RXS spectrum at the first satellite spot 
consists of the pure rank one energy profile (a^-j'(LL')) in the uniform helical phase. Upon 
cooling temperature, the spectrum comprises that of the rank two {a^^^{Lo)) as well as a^l{uj) 
both in the spin slip and conical phases. On the other hand, the spectrum at the second 
satellite spot is made up of the pure a^-[{uj) in the uniform helical and conical phases. In 
the spin slip phase, however, the spectrum includes a^l{Lo) too. Our numerical calculation 
have revealed that the magnitude of the |a£-j'(u;)| is much larger than for instance, 

|ag^^(LL')|/|agj'(a;)| ~ 27.7 with T = 0.6 eV. Thus, whenever both profiles mix, we can expect 
the a^^{uj) dominates the spectral shape. Then, the spectrum at the first satellite spot looks 
like |a£-j'(u;)p at all magnetic phases. Contrary to this, we conclude that the spectral shape 
of the second satellite spot varies as a function of temperature. That is, the spectrum is pure 
|a!g-i^ (aj)p in the helical phase, then looks like jaj^^-j' in the spin slip phase, and finally, 
becomes pure |Q;g^^(a;)p again in the conical phase. If such variation of the spectral shape 
would be observed by the experiment, it demonstrates a clear evidence of the change of the 
magnetic structures. 

If the experiment in the conical phase is attainable, another attractive outcome can be 
expected. That is, the RXS spectrum from the higher-order harmonic satellite spots may be 
detected. It is brought about by the distortion of the orientation of the magnetic moment 
within the basal plane. In neutron scattering, the signals at the (6 it l)th magnetic satellite 
spots were observed. ""^ In the present case of RXS in the El transition, the intensities at the 
(6 lb 2)th as well as (6 ib l)th order satellite spots, which are related to rank two and one 
operators, respectively, will be detected. By evaluating the Ho case in the vicinity of the M5 
absorption edge, we can conclude that the intensity at the fourth satellite spot is within the 
reach of the present experimental condition of the El transition in the tt — a' channel. 

As for the numerical results, first, we have concentrated on the comparison of our result 
with the experiment in the uniform helical phase. The calculated RXS spectrum at the first 
satellite spot and the absorption coefficient show excellent agreement with the observed ones 
both the spectral shape and the peak position (for the former spectrum).^' ^^'^^ The agreement 
assures the reliability of the present analysis based on the localized picture. At the same time, 
it also demonstrates the effectiveness of the spectral shape analysis in the RXS theory, which 
is rare, in particular, in the /-electron systems. The agreement of the spectrum at the second 
satellite spot between ours and the experiment is not an ideal one. Although the peak position, 
below about 1.8 eV compared to that of the first satellite spectrum, is properly reproduced, 
the theory includes an extra peak around 1353.2 eV which is missing in the experiment. The 
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reason of this discrepancy is still unclear. 

Next, our evaluation of the spectra at the M4 edge showed that the magnitudes of them 
are roughly two orders of magnitude smaller than those at the M5 edge. Thus, the possibility 
of experimental detection seems to be slim at the second satellite spot. On the other hand, 
the magnitude of the spectrum at the first satellite spot at the M4 is expected to be the same 
order as that of the second satellite spot at the M5 edge, and seems to be detectable. Finally, 
we have confirmed that the spectrum the second satellite spot in the spin slip phase certainly 
governed by the |Q!g^^(u;)| as anticipated. 

A qualitative part of the present analysis may be applicable to the interpretation of the 
RMXS experiment near the Ho L3 edge in the El process, by neglecting the band nature of 
the 5d electrons as the previous works had followed.^'^'^"^'"^^ Since the absorption coefficient 
at the L3 edge exhibited no multi-peak structure, we conclude either the multiplet splitting 
of the intermediate are not large enough and/or the energy resolution were not fine enough 
to distinguish the line shapes of the different profiles a^l{uj) and a^^^uj). Thus, only we 
can guess is the polarization analysis. If the sample used were in the uniform helical phase, 
eqs. (6.3) and (6.4) state that the spectra at the first and second satellite spots are pure 
a^l{uj) and a^^{uj), respectively. The former is absent in the a — a' channel. The latter has 
contributions both in the a — a' and a — n' channels. These conclusions are just the ones 
deduced previously."^^'"^^'"^^'^^ Note that, in ref. 13, although the authors concluded there was 
no El contribution in the a — tt' channel at the second satellite spot (0, 0, 2 + 2r), the data 
exhibits a small hump around 8071 eV, at which the El peak is expected. By substituting 
the value of the Bragg angle (6^'^^ ~ 18°) in eq. (6.4), the intensity of the a — tt' channel in 
the El process is evaluated about a tenth of that in the a — a' channel. Since the ratio is 
nearly equal to that observed, the data may show the contribution from the El process in the 
a — it' channel at the second satellite spot. Thus the intensity around this energy may show the 
resonant behavior. Additionally, the analysis on the E2 process by means of the present theory 
qualitatively shows complete agreement with the previous interpretation.^'^' ^^'^^ Because our 
formalism allows naturally inclusion of the energy profile analysis, such an investigation is the 
next step along this line.^^ 

Since the wavelength of soft x-ray is suitable for the investigation of the long periodic 
magnetic structure observed in several heavy rare earth material, the present theory may be 
effective for the analysis of the R(M)XS spectra at the M^^^ edges of such heavy rare earth 
compounds, for example, DyB2C2. Actually, the absorption coefficient in the vicinity of the 
Dy M5 edge exhibits multi-peak structure similar to that observed in Ho.^'^^ On the other 
hand, the RXS spectra differ significantly between those of Dy and Ho. The spectrum of Dy 
extends over a few tens eV, an order of magnitude wider than that of Ho, which may not be 
explained within the present localized electron picture. The work along this line is a future 
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Appendix: The geometrical factors and the scattering amplitude 

Here, we present the geometrical factors needed to evaluate cqs. (5.10) and (5.11). The 
results are classified by the photon polarizations. The origin of ip is defined such that the y 
axis lies in the scattering plane. The geometrical configuration we adopted in this work can 
be found in our previous paper. For the first satellite spot Gi = (0, 0, ri). 



(e' X e)^-z(e' x e)^, = e^^cos^^ < 
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where 0^^^ is the Bragg angle associated with Gi. Then, for the second satellite spot G2 
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where 9^'^^ is the Bragg angle associated with G2. Note that a relation sin^(^) = 2sin^(-'^) 
holds. 

By substituting these results into eqs. (5.10) and (5.11), we obtain the RXS intensities at 
the first satellite spot as 
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These results show the spectral intensities exhibit no azimuthal angle dependence. 
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